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1 Abstract 

A generic proecological traffic control model for the urban street canyon is proposed by de- 
velopment of advanced continuum field hydrodynamical control model of the street canyon. 
The model of optimal control of street canyon dynamics is also investigated. The mathemat- 
ical physics' approach (Eulerian approach) to vehicular movement, to pollutants' emission, 
and to pollutants' dynamics is used. The rigorous mathematical model is presented, us- 
ing hydrodynamical theory for both air constituents and vehicles, including many types of 
vehicles and many types of pollutants emitted from vehicles. The six proposed optimal mon- 
ocriterial control problems consist of minimization of functionals of the total travel time, of 
global emissions of pollutants, and of global concentrations of pollutants, both in the studied 
street canyon, and in its two nearest neighbour substitute canyons, respectively. The six op- 
timization problems are solved numerically. Generic traffic control issues are inferred. The 
discretization method, comparison with experiment, mathematical issues, and programming 
issues are discussed. 

2 Description of the model. 

In the present article we develop a continuum field model of the street canyon. In the 
next article we will deal with numerical examples [1]. The vehicular fiow in the canyon 
is multilane bidirectional one-level rectilinear, and it is considered with two coordinated 
signalized junctions [2, 3, 4]. The vehicles belong to different vehicular classes: passenger 
cars, and trucks. Emissions from the vehicles are based on technical measurements and 
many types of pollutants are considered (carbon monoxide CO, hydrocarbons HC, nitrogen 
oxides NOx). The vehicular dynamics is based on a hydrodynamical approach [5]. The 
governing equations are the continuity equation for the number of vehicles, and Greenshields' 
equilibrium speed-density u-k model [6]. 
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The model of dynamics of pollutants is also hydrodynamical. The model consists of a set 
of mutually interconnected nonlinear, three-dimensional, time-dependent, partial differential 
equations with nonzero right-hand sides (sources), and of boundary and of initial problem. 
The pollutants, oxygen, and the remaining gaseous constituents of air, are treated as mixture 
of noninteracting, Newtonian, viscous fluid (perfect or ideal gases). The complete model 
incorporates as variables the following fields: density of the mixture, mass concentrations of 
constituents of the mixture, velocity of mixture, temperature of mixture, pressure of mixture, 
intrinsic (internal) energy of mixture, densities of vehicles, and velocities of vehicles. The 
model is based on the assumption of local laws of balance (conservation) of: mass of the 
mixture, masses of its constituents, momentum and energy of the mixture, the numbers of 
the vehicles, as well as of the state equations (Clapeyron's law and Greenshields' model). 
The equations of dynamics are solved by the finite difference scheme. 

The six separate monocriterial optimization problems are formulated by defining the 
functionals of total travel time, of global emissions of pollutants, and of global concentrations 
of pollutants, both in the studied street canyon and in its two nearest neighbour substitute 
canyons. The vector of control is a five-tuple composed of two cycle times, two green times, 
and one offset time between the traffic lights. The optimal control problem consists of 
minimization of the six functionals over the admissible control domain. 

3 Equations of dynamics. 

Under the above model specifications, the complete set of equations of dynamics of the model 
is formulated as follows (we follow the general idea presented in [7, 8]): 

El. Balance of momentum of mixture - Navier Stokes equation. 



293.16 [K]), ^ is the second viscosity coefficient (^ = 15.6 ■ 10~^[^] for air at temperature 
T — 293.16[K]), F = pg is the gravitational body force density, g is the gravitational 
acceleration of Earth (g = (0, 0, — 9.81)[^]), Vv is gradient of the vector (so it is a tensor of 
rank 2). We assume that the gaseous mixture is a compressible and viscous fluid. 
E2. Balance of mass of mixture - Equation of continuity. 




(1) 



dp 
di 



+ div(pv) = S. 



(2) 



We have assumed the source Dl. 

E3. Balances of masses of constituents of mixture - Diffusion equations. 



E3a. 



(3) 
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N-1 u 
+ E {(Am - Aat) • div[pV(c^ + ^VT)]}, 1=1,..., Ne. 



m=l 



E3b. 



p(^ + voVq) = -q5+ (4) 
+ 5: {(Am - Aiv) • div[pV(c^ + ^VT)]},i = (7V^ + 1), TV, 

m=l 

where A^ = -Dmi is the mutual diffusivity coefficient from the i-th constituent to m-th one, 
and Ai is the autodiffusivity coefficient of the z-th constituent, and /cT,m is the thermodiffu- 
sion ratio of the m-th constituent. The diffusivity coefficients and thermodiffusion ratios are 
constant and known). In E3a we have assumed the sources D1-D2. In E3b only the source 
Dl is taken into account. Since the mixture is in motion, wc cannot neglect the convection 
term: v o Vq. We assume that the barodiffusion and gravito diffusion coefficients are equal 
to zero. 

E4. Balance of energy of mixture. 

de 1 
PijT. + V o Ve) = -(--v^ + e),5 + T : Vv + div(-q) + a, (5) 

where e is the mass density of intrinsic (internal) energy of the air mixture, T is the stress 
tensor, symbol : denotes the contraction operation, q is the vector of flux of heat. We assume 
that [2]: 

N 

e, = J-{c.fc^Texp(-^) • [{-'-) • (1 - exp(-^)) - exp(-^)]} + (7) 
rrii kb-l c Kb-L KbJ- 

^ - ^, (8) 

m 

/X, = kBT . {ln[{c,p){kBT)-^{^){^)^} + m,\g\z, (9) 

Tmfe = -pSmk + (10) 

' ^ ~ \^"'''^^^^^''''dx^'^ ^ ^ ' l.i^"^kdiv{v)f],m, k = l, 3, 

T:Vv=EETm.|^, (11) 

m=l k=l ^-^k 

^ BT 

q = E{[(7^ + /i^)J^] + Mvr]}, (12) 
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j, = -pAi(c. + ^VT), (13) 
a,, = ''P^ , (14) 

A = [pA.] • - .ly^"^-"-""^ J , (15) 

where is the mass density of intrinsic (internal) energy of the ith constituent of the air 
mixture, rrii the molecular mass of the 2th constituent, ks = 1.3807- 10^^'^fj^] is Boltzmann's 
constant, /Xj is the complete partial chemical potential of the ith constituent of the air mixture 
(it is complete since it is composed of chemical potential without external force field and of 
external potential), mgir = 28.966 [u] is the molecular mass of air (l[u] = 1.66054-10"^'' [kg]), 
is Kronecker's delta, Cp^i is the specific heat at constant pressure of the zth constituent of 
air mixture, h = 6.62608 - 10"^"^ [J - s] is Planck's constant, jj is the vector of flux of mass of 
the ith constituent of the air mixture, and k, is the coefficient of thermal conductivity of air. 
These magnitudes are derived from Grand Canonical ensemble with external gravitational 
Newtonian field. 

E5. Equation of state of the mixture - Constitutive equation - Clapeyron's 
equation. 

P^^.T (16) 

is Clapeyron's equation of state for a gaseous mixture, where R — 8.3145 [ ^^^^ -^ ] is the gas 
constant. 

Pi^Ci p (17) 

m 

are partial pressures of constituents according to Dalton's law. 

E6. Balances of numbers of vehicles - Equations of continuity of vehicles. 

-^+dw{kl,wlJ^O. (18) 
E7. Equations of state of vehicles - Greenshields model. 

wf,.,(x, t) = (wl,j - (1 - %^), 0, 0). (19) 

The Greenshields equilibrium speed-density u-k model is assumed [6]. The values of 
maximum free flow speed wf^jj, and of jam vehicular densities kf^^.^^^, are given in Tables 
3 and 8 of [1]. 

E8. Technical parameters. 

The dependence of emissivity on the density and velocity of vehicles is assumed in the 
form: 
E8a. 
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^i,ct,vt\^^ ^) = Kvt\^^ t) ■ [(-7^ ) • {ect,vt,ii+i - ect,vt,H) + ect,vt,iiU (20) 

where Wct,vt,H are experimental velocities, \y^f^^^{x,t)\ G {wct,vt,ii^'^ct,vtdi+i)^ ect,vt,ii, are ex- 
perimental emissions of the ctth exhaust gas from single vehicle of vtth type at velocities 
Wct,vt,iii respectively, measured in [:;;^], ii — 1, Nemi Nem is the number of experimental 
measurements. Similarly, the dependence of the change of the linear density of energy on 
the density and velocity of vehicles is taken in the form: 
E8b. 

(^l,vA^^ t) = Qvt ■ \vti^, t) ■ [(— T -T ) • [(T^t,H+l - (^vIm) + (^vtMl (21) 

where <T„t,ip are experimental values of consumption of gasoline/diesel for a single vehicle of 
vtth type at velocities tw^jt^j,, respectively, measured in [;;;^], Qvt is the emitted combustion 
energy per unit mass of gasohne/diesel [j^]. 



4 Optimization problems. 

Our control task is the minimization of the measures of the total travel time (TTT) [5], 
emissions (E), and concentrations (C) of exhaust gases in the street canyon, therefore the 
appropriate optimization problems may be formulated as follows [2]: 
FO. Vector of control. 

u=(<7i,Ci,^2,C2,F)e[/^'i-, (1) 

where u is vector of boundary control, gm are green times. Cm are cycle times, F is offset 
time, and f/^<^™ is a set of admissible control variables (compare A8, A9, B5, B5S, B5SS). 

We define six functionals F1-F6 of the total travel time, emissions, and concentrations 
of pollutants in single canyon, and in canyon with the nearest neighbour substitute canyons, 
respectively. 

Fl. Total travel time for a single canyon. 

Jttt(u) = EEE/ L KvtMdxdt. (2) 
F2. Global emission for a single canyon. 

2 Us CT VT Ts 

^e(u) = EEEE / <ct,.t(^, t)dx dt. (3) 

s=l 1=1 ct=liA=l-'^ •'^ 
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F3. Global pollutants concentration for a single canyon. 

^B^l ra i-b i-c i-Ts 

Jc(u)=psTP- J2 Ci(x,y,z,t)dxdydzdt. (4) 

■^-^^ Jo Jo Jo Jo 

F4. Total travel time for the canyon in street subnetwork. 

^TTT,ext(u) = JtTt(u) + (5) 

2 Us VT 

+a ■ J2 "TTT,ext E E ^i%*,jam " (^s " ^.)- 
s=l 1=1 vt=l 

F5. Global emission for the canyon in street subnetwork. 

</E,ext(u) = Je(u) + (6) 
2 Us CT VT 

+a ■ ^ a^^g^t JL^Yl ef,ct,^t,ja™ • (Cs - 9s)- 

s=l 1=1 ct=l vt=l 

F6. Global pollutants concentration for the canyon in street subnetwork. 

^C,ext(u) = Jc(u) + (7) 

Ne-1 2 

+PSTP -a-b-C- Ci.STP • ac,ext • (Cs - gs)- 

i=l s=l 

The integrands ki^vfi^i,ct,vf^^ functionals F1-F6 depend on the control vector u FO 
through the boundary conditions B0-B8, through the equations of dynamics E1-E8, as 
well as, through the sources D0-D2. The value of the vector of control u directly affects 
the boundary conditions B5, B5S, B5SS, and then the boundary conditions B6-B8 for 
vehicular densities, velocities, and emissivities. It also affects the sources D0-D2. Next, 
it propagates to the equations of dynamics E1-E8 and then it influences the values of 
functionals F1-F6. We only deal with six monocriterial optimization problems 01-06, 
and not with one multicriterial problem. We put the scaling parameters equal to unity: 
'^TTTcxt ~ '^Ecxt ~ "^Ccxt ~ ^-^i fuuctiouals F4-F6. psTP is the density of air at standard 
temperature and pressure STP, q^stp is concentration of the ith constituent of air at standard 
temperature and pressure. Jttt and JxTT.ext are measured in [veh • s], Je and JE,ext are 
measured in [kg], and Jc and Jc,ext are measured in [kg ■ s], respectively. 

Now we formulate six separate monocriterial optimization problems 01-06 that consist 
in minimization of functionals F1-F6 with respect to control vector FO over admissible 
domain, while the equations of dynamics E1-E8 are fulfilled. 

01. Minimization of total travel time for a single canyon. 

J*^:p = Jttt(uttt) = minju e U""^^ : Jttt(u)}; (8) 

02. Minimization of global emission for a single canyon. 

^E = ^e(u^) = min{u e U^"-^ : Je(u)}; (9) 
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03. Minimization of global pollutants concentration for a single canyon. 

= Jc«) = min{u e U'^^ : Jc(u)}; (10) 

04. Minimization of total travel time for a canyon in street subnetwork. 

^TT.ext = -^TTT.extlu^TT.ext) = min{u G t/"^"' : JTTT,ext(u)}; (11) 

05. Minimization of global emission for a canyon in street subnetwork. 

4ext = ^E,ext(u^,ext) = niin{u G t/^^^ : JE,ext(u)}; (12) 

06. Minimization of global pollutants concentration for a canyon in street 
subnetwork. 

^^,ext = ^c,ext(u*c,ext) = min{u G U^"^ : Jc,ext(u)}, (13) 

where J^tt? ^e^ Jqi "^TTTcxti ^Ecxt^ ^Ccxt^ ^he minimal values of the functionals F1-F6, 
and u^TT) ^E) '^C' ^TTT,ext) ^E,ext) ^c.ext' control vectors at which the functionals reach 
the minima, respectively. 

5 Computer simulations. 

Prom thousands of performed optimizations 01-06 wc selected some rcprcscntable opti- 
mizations. The C language programme of more than 14000 lines written by the author was 
being run on workstations of Hewlett-Packard, Sun, Silicon Graphics, under UNIX operating 
system. It was being run also on Pentium personal computers under Linux operating sys- 
tem. We used C language UNIX/Linux compilers cc. The time of simulations varied from 
a couple of hours to several days or even weeks, depending on the discretization parameters. 
The discretization steps were tested in order to obtain finite solutions to the optimization 
problems in reasonable time of simulation. They were also studied from the point of view of 
sensitiveness of optima appearing on them. 

6 Mathematical questions. 

The problem of uniqueness of solutions of optimization problems 01-06 is complex. If 

there are no vehicles allowed in the canyon, then all the six functionals 01-06 are constant 
functions of parameters, the optima exist, and the number of optimal solutions is continuum. 
The values of functionals are then zero for total travel time and emissions, and are positive 
for concentrations. These results are analytical. Por other cases there are no such analytical 
results for both the uniqueness and existence known to the author. The numerical solutions 
of optimization problems are approximately global within the error induced by discretization 
of the physical domain S and of the domain of control parameters U^*^^, The method of 
optimization was a full search. 
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7 Conclusions. 



The proecological traffic control idea and advanced model of the street canyon have been 
developed. It has been found that the proposed model represents the main features of 
complex air pollution phenomena. 
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